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ABSTRACT 

Energy spreading loss (ESL) is qualitatively defined as the reduction in peak 
echo level due to energy spreading of the transmitted acoustic pulse in time. An 
analysis of the impact of shallow water propagation on ESL was performed with the 
aid of a high performance computer using the FEPE_S YN and EXT_TD programs to 
compute the spreading of the received pulse due to multipath propagation in shallow 
water. A Blackman windowed pulse was used to model the transmitted pulse, which 
was centered at 3.5 kHz, with 200 Hz bandwidth. For input parameters, typical 
seasonal sound speed profiles and a Hamilton geoacoustic model of Area Foxtrot off 
the U.S. eastern seaboard was used. ESL’s impact on sonar performance was 
determined as a function of range, source and target depth, sound speed profiles and 
geoacoustic properties. The impact of shallow water propagation on the correlation 
of the transmitted and propagated pulses through the quantitative definition of 
mismatch loss (MML) was also discussed. 

The results showed that strong ESL (5 to 10 dB) existed over a sand (reflective) 
bottom and was generally invariant with range. ESL was correlated with TL, i.e., 
areas of high spreading loss were found in regions of high TL. ESL was not as large 
(3 to 5 dB) over silt/clay (absorptive) bottoms due to the increased absorption of the 
bottom refracted path thus reducing the number of multipath modes. Broadband 
pulses were found to exhibit fewer fluctuation than single fi'equency signals, and 
generally the total TL loss was a few dB larger than a single cw case. To overcome 
the ESL, integration techniques based on an accurate prediction model in the post 
analyzing system are required with a high temporal resolution of the echo energy 
shape. 
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I . INTRODUCTION 

A. NAVY'S INTEREST IN ASW 

After the end of the cold-war the U.S. Navy's strategic 
interests shifted from deep water to shallow water issues 
because of the decrease of the threat of global war and the 
increased threat of regional conflicts. It is anticipated that 
these smaller scale conflicts will probably occur in littoral 
regions, not in the deep ocean. Since there is no longer a 
large foreign Navy to challenge the U.S. Navy on the open 
seas, regional conflicts, dominated by quiet diesel submarines 
and densely populated minefields, have become the USW 
challenge of the present. 

The seas of littoral regions are usually shallow and 
their high spatial and temporal variability requires a high 
degree of knowledge of shallow water oceanography and 
acoustics to counter the ASW and mine warfare threats. Shallow 
water regions by definition are very complex oceanographically 
and acoustically due to temporal, spatial and spectral 
variations of environmental parameters. 

Recent developments in sound silencing technology have 
resulted in submarines becoming very quiet, and advances in 
air injection propulsion (AIP) systems have made very long 
operation times on the battery possible. These improvements 
have increased the importance of active sonar operations in 
shallow waters. Accordingly, the assessment and improvement of 
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active sonar performance in shallow water is now a major 
concern to planners and operators . 

B . BACKGROUND 

Shallow water acoustics has long been a topic of ASW 
research but with emphasis on low frequency (< 1000 Hz) 
propagation [Urick, 1983] . Observations and model studies have 
shown that low frequency sound propagation in shallow water, 
especially for situations where the geoacoustic impact of the 
bottom/ subbottom is significant [Hamilton, 1972, 1979, 1980, 
1985, 1987; Jackson and Briggs, 1992; Hamilton and Bachman et 
al, 1982; Mourad and Jackson, 1983], is a complex problem. 
The complexity arises because many environmental acoustic 
parameters, such as the tenperature and salinity of the water 
column, bathymetry, type of sediment, background noise, marine 
life, etc.., are highly variable, both temporally and 
spatially. 

For tactical active sonar performance, the focus is on 
higher frequency sound propagation (order of several kHz) , 
which is generally addressed in terms of ray acoustics. Ray 
acoustics is an accurate approximation to the full wave 
equation at high frequencies in deep water where diffraction 
effects are negligible. In shallow water diffraction and 
multipath interaction limit the use of ray acoustics and full 
wave solutions are required [Scanlon, 1995] . Although a full 
wave equation solution requires a significant computer 
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intensive effort at high frequencies, accurate approximations 
such as the Finite Element Parabolic Equation (FEPE) [Collins, 
1988] exist to solve this problem and it is this model which 
has been selected for use in this analysis. FEPE is a single 
frequency (i.e., zero bandwidth) model for a continuous wave 
source. A time domain version of FEPE, called FEPE_SYN 
[Collins, 1989] , is necessary to predict the time arrival 
structure of a pulse of finite bandwidth. 

Because of the small range and depth steps required at 
the relatively high frequency of 3.5 kHz, the calculation of 
the full wave solution would have taken an astronomical length 
of time several years ago. Now computing power is available 
to realistically process time series data with modern 
powerful computers . 

It is difficult to improve sonar performance in shallow 
water because of the phenomena of “Energy Time Spreading 
This phenomena has been known for more than 20 years [Bell, 
1990] as it relates to low frequency passive acoustics, and 
was reported as time stretching due to “multipath effects". 
However, energy spreading loss (ESL) has not been widely 
modeled or measured as an active sonar equation parameter. 
Jensen (1993) and several NUSC researchers [Jones, 1990; Bell, 
1990] have performed some excellent initial research on ESL 
for active acoustic sonar systems recently. Since their 
research efforts were based on observational data, no ESL 
modeling capability was available which attempted to solve the 
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full wave equation in the time domain. 

In this research, ESL is treated as an active sonar 
equation parameter as an additional loss in addition to the 
standard terms. The models “FEPE_SYN" [Collins, 1989] and 
“EXT_TD" [Rovero, 1992] are used to estimate its magnitude and 
time spreading characteristics based upon an acoustic pulse of 
finite bandwidth transmitted from an active sonar system in 
shallow water. EXT_TD is used to add signal pulses of finite 
bandwidth to the time domain model FEPE_SYN, 

1. Definition of ESL 

Energy spreading loss (ESL) is defined qualitatively as 
the reduction in peak energy (or power) level due to the 
spreading of the energy of a transmitted acoustic pulse with 
time beyond the original transmitted pulse length. ESL is 
quantitatively defined in two ways: the first is the reduction 
in peak energy level to total energy [Jones, 1990] ; the second 
is the coherence loss of a matched filter due to mismatching 
of the pulse shape [Jensen and Sabbadini, 1993] . The first 
energy- based definition of ESL is termed “ESL” while the second 
coherence measure is termed the "mismatch loss (MML)" to 
differentiate them physically. 

ESL is based on the time stretching of the transmitted 
pulse, in which the pulse is stretched in time by multipath 
propagation. The time stretching causes the peak energy of the 
received pulse at a given point to be reduced below that of an 
echo at the same range but for which no multipath propagation 
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occurs, e.g., in deep water. 

MML is based on the change in correlation or coherence 
between the transmitted and received pulse shapes at the 
receiver, a definition proposed by Jensen and Sabbadini (1993) 
for the deep bottom reflected case. 

2. Past ESL Results 

ESL has been studied over the past 20 years by a number 
of researchers as a multipath effect [Bell, 1990] . ESL is a 
complex phenomena of multipath or multimode wave propagation 
prevalent in shallow water. The term “energy splitting loss” 
was introduced by Stewart and Brandon (1967) to describe the 
type of signal distortion process whereby multipaths can split 
the echo energy into a number of resolvable arrivals. 

“Energy spreading loss” was suggested by Weston (1965) 
who observed that the correlation loss due to spreading in the 
time delay was precisely the same phenomenon as the time 
spread and associated reduction in peak level observed when a 
short ping was transmitted. 

The standard deviation, o, of the stretched pulses in the 

time domain about the mean arrival time of the pulses has been 
used as a measure of time stretching and has been widely used 
in previous ESL studies [Van Trees, 1971] . The following 
equations summarize the use of the standard deviation as a 
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parameter to measure ESL: 



A^(t)dt = 51 [s * 



spread (k)] 



( 1 - 1 ) 



o 



= [ J (t-fi)' A'(t)dt 

= 51 [ (k*dt - * S * A^pread(k) ] 



( 1 - 2 ) 



where fj, = 




spread (k)] . 



( 1 - 3 ) 



Here A^ represents normalized power, t represents time [sec] , 



k th time increment of the stretched pulse, k is the index 
number of the time increment and s is the normalization noise 
factor. Unfortunately, this measure of ESL was found to be 

corrupted by ambient noise in measurements of o during at-sea 
experiments . 

Bell ( 1990 ) modified Weston's definition of ESL by 
incorporating a Rayleigh distribution instead of a Gaussian 
distribution for the spread pulses. He used a Monte Carlo 
technique to include the effects of Rayleigh fluctuations in 
each resolvable part of the signal return. Bell's work was 
extended by Jones ( 1990 ) and his empirical equation (termed 
the "Bell -Jones Equation") gave good agreement with observed 
data [Chan, 1992 ] . Because the Bell-Jones equation related the 



H, the mean arrival time, A^gpread(k), the acoustic power at the 
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ratio of a time spreading function and system resolution, 
o /R, to ESL, a threshold was set above a specified ambient 

noise level and a Gaussian distribution of the time spreading 
was assumed. Jones related the ratio of total energy (energy 
computed over the entire stretched signal) to peak energy 
(peak signal power times the resolution cell width) to ESL 
[Young, 1988] , 

ESL = 10*logio (TE / PE) [dB] . (1-4) 

where TE represents the total energy and PE the peak energy. 
A resolution cell width, R, typically would be defined as the 
reciprocal of the correlator bandwidth BW [Hz] . 

R = 1 / BW [sec] . (1-5) 

The ESL defined above is termed the “Jones ESL" in this study, 
which is used as a second measure of ESL. 

Recently, Chan (1992) concluded that the multipath 
boundary interaction in the shallow water propagation channel 
accounted for most of the loss in expected performance of the 
wideband LPM (coded pulse) waveforms implemented on the 
AN/SQC-53C sonarset . The results indicated that ESL was a 
major loss mechanism in shallow water, and up to 12 dB of 
degradation was measured for specific waveforms. A modeling 
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technique which incorporated all frequencies in the signal 
band was used to model the time spreading in this experiment. 
Unfortunately, the model used is not considered to be accurate 
in shallow water. 

Jensen and Sabbadini (1993) proposed MML (mismatch loss) 
as a measure of the signal degradation experienced by low 
frequency active (LFA) sonars. MML is based on the change of 
the pealc pulse shape at the receiver compared to the shape of 
the transmitted pulse. The normalized coherence between the 
replica (transmitted) and a propagated pulse is given by, 



where Vr(t) is the replica pulse, Vp(t) is the propagated or 
received pulse near the maximum peak of the received pulse, 
and the interval for integration is usually taken as the 
duration of the replica pulse for matched filter processing 
[Jensen and Sabbadini, 1993] . MML is considered the third 
quantitative measure of ESL. Jensen and Sabbadini 
investigated LFA bottom bounce losses using a bottom 
interaction simulator and 1/100 scale cylinder -shaped target 
with rounded end caps. The bottom impulse response functions 
were obtained from SACLANT CENTER data. The experiment was 




(1-6) 
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based on a bistatic geometry, and two bottom bounces for a 
pulse emitted at various launch angles, including 0, 45, 90 
degrees of target aspect. Their results indicate that MML can 

reach 5.3 dB for a wideband (0.5 ~ 5 kHz), 1.5 sec LFM pulse 
(without target) in matched filter processing. Only under 
ideal (no spreading) conditions will LFA sonar performance 
predictions be realized when coherent processing is used. 

In this study ESL and MML are calculated for one way 
propagation only (without the target) to focus the research on 
the impact of the environmental acoustic parameters affecting 
propagation on ESL/MML. A windowed Blackman pulse is used in 
all model runs to represent the outgoing pulse shape. This 
study will not address the impact of the transmitted wave 
pulse type, duration, etc., on ESL and MML. 

The model results from this research are consistent with 
the measured ESL results of Chan (1992) , where time spreading 
and energy spreading were the product of two way (monostatic) 
TL plus target scattering loss. It was not possible for Chan 
to analyze many pulses with exactly the same propagation 
environment since measured data (containing echo and 
reverberation) were used in his analysis, and thus he was 
forced to average different phased signals statistically. 
Similarly, Jensen and Sabbadini's analysis was limited because 
it did not account for multipath propagation, and therefore 
their MML was based on only one ray or predominant mode of the 
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acoustic wave. 



3 . Impact of ESL on Tactical Active Sonar Performance 

Although ESL is generally incorporated as an additional 
loss term added in an ad hoc fashion to the active sonar 
equation, it is predicted by the wave equation. Without 
recognizing it, ESL has been measured as a part of TL in many 
past measurements [Urick, 1903] and is a major factor in total 
energy loss to the beamformer (excluding geometrical 
spreading) in shallow water. Whenever the wave equation is 
solved exactly, the solution includes ESL as well as TL if the 
model results are defined properly. Since the FEPE model 
solves the wave equation accurately, FEPE estimates ESL 
accurately. In this study a range independent environment is 
assumed that is invariant in time and the significant problems 
of sound speed fluctuations, sloping bottoms and spatially 
variable geoacoustic properties in shallow water are avoided. 

For tactical active sonars, ESL can exert a significant 
degradation on sonar performance by reducing the detection 
range significantly in shallow water. This research 
demonstrates the impact of ESL degradation on active sonar 
performance as a function of source and target depth, range, 
sound speed profile and geoacoustic properties of the bottom. 

Clearly TL is due to a geometrical spreading loss, 
absorption, and other boundary losses prior to the energy 
reaching the array. This energy can never be regained at the 
receiver location. However, it is possible to partially regain 
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ESL since the ESL energy arrives at the target location but is 
spread in time well beyond the original pulse width. This 
means it is possible to improve sonar performance if we can 
accurately model the spread of the transmitted pulses, and 
apply advanced beamforming/signal processing techniques such 
as Inverse Beamforming (IBF) [Nuttall and Wilson, 1991; 
Wilson, 1995; Fabre and Wilson, 1995]. 

In deep water the bottom boundary effects on propagation 
are relatively small and ESL is usually negligible at short 
range (< 15000 m) . In shallow water where the bottom boundary 
effects are significant, ESL is large. Multipaths or 
multimodes generated by boundary interactions cause 
significant time stretching and result in large ESL. The 
physical description of time spreading which degrades sonar 
performance in shallow water is discussed in the next chapter. 



C . OBJECTIVE 

This study attempts to define ESL quantitatively as a 
function of several parameters, e.g., range, source and target 
depth, water sound speed profile and composition of the sea 
floor. An analysis of the impact of shallow water propagation 
on ESL (based on this definition) is performed. The influence 
of ESL on sonar performance is discussed and recommendations 
are made to overcome the degradation due to ESL. An 
examination of the previous ESL definition [Jones, 1990] is 
also performed. 
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A second, but minor, objective is to analyze the impact 
of shallow water propagation on the correlation of the 
transmitted and propagated pulses through the quantitative 
definition of MML [Jensen and Sabbadini, 1993] . 



12 



II. MODELING ESL IN SHALLOW WATER 



A. SIGNAL PROPAGATION MODELS USED IN THIS ANALYSIS 

1. FEPE_SYN and EXT_TD 

FEPE_SYN is a fully range - dependent , time domain version 
of FEPE [Collins, 1989] in which the ocean transfer function 
(similar to TL) is calculated as the output. FEPE_SYN is 
implemented with a "window” in range and depth to limit the 
area over which the output is produced and a logfile is 
created which provides the structure of the data. FEPE_SYN was 
augumented by Rovero (1992) with a software program called 
EXT_TD to incorporate the propagation of finite length pulses 
with broadband responses and can easily be expanded to 
accommodate any new parameters required to match the 
transmitted pulse characteristics of existing tactical active 
sonars . 

EXT_TD is a separate program which reads the output of 
FEPE_SYN and creates a time domain signal at a selected range 
and depth. A minor modification was made by the author to 
compare the resultant time spreading of the output pulse to 
the transmitted pulse length in order to calculate ESL values 
defined in the next section. As stated earlier, the source 
signal pulse shape used in this analysis is a Blackman 
windowed pulse. 

Pulse propagation modeling is concerned with simulating 
the effects associated with the transmission of a finite 
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bandwidth signal characterized by a known frequency spectrum 
as contrasted with the more familiar single f reqfuency 
continuous wave propagation. In principle, the frequency 
domain wave equation treats broadband signals by Fourier 
synthesis of individual cw solutions over the frequency 
spectrum. In the presence of nonlinearities, however, 
interaction among frequency components invalidates this 
frequency domain approach. In the time domain, the wave 
equation can be formulated using methods which remove 
pathological limitations from numerical solutions. [Etter, 
1991] 

Computers used in this research are the Power Onyx (2.75 
MHz MIPS R8000CPUS, 64 bit operating system) Silicon Graphics 
Computer and the Model J916/4 Cray Computer, a powerful Cray 
research cortputer. It is impractical to obtain the ESL results 
without these, or even more powerful, computers at 3.5 kHz 
using a full wave propagation model such as FEPE_SYN. 

2 . Oceanographlc/Geoacoustlc Inputs 

The input parameters were selected for a typical tactical 
sonar having a center frequency of 3.5 kHz, bandwidth of 200 
Hz and at a depth of 7.3 m. An example of a typical FEPE_SYN 
input file is shown in Table 1. Oceanographic inputs were 
seasonal sound speed profiles and bathymetry profiles of Area 
Foxtrot, a tactical active sonar exercise area, south of Long 
Island. The geoacoustic inputs were a full Hamilton 
geoacoustic representation of sand and silt/clay sediment 



14 



types [Scanlon, 1995] found in Area Foxtrot. 

3 . FEPE 

FEPE calculates the transmission loss only for a single 
frequency [Collins, 1989] and the latest version of FEPE was 
used to illustrate the TL as a function of range and depth in 
this study. Because ESL is strongly associated with TL (shown 
in Chapter V) and must be calculated separately from it, the 
FEPE TL results for a single frequency are used for comparison 
of TL from FEPE_SYN for the full 200 Hz signal bandwidth. 

4. Strategy for Varying Input Parameters 

For a frequency of 3.5 kHz, FEPE_SYN is very computer 
intensive. Sait^^ling theory requires that a very small temporal 
increment be selected to avoid aliasing. The finite element 
method also requires a very small range and even smaller depth 
mesh size, to avoid spatial aliasing. In order to generate 
acoustic pulses with a maximum frequency of 3600 Hz, the 
Nyquist frequency of 7200 Hz is a minimum frequency limit for 
analysis. To be conservative, we chose a sampling rate of 4 * 
3600 Hz = 16384 Hz to avoid aliasing. A 16,384 point FFT was 
used, producing a 1 Hz frequency resolution for a 1 sec time 
window. 

The smallest wavelength in the transmitted pulse was used 
to determine the range/depth mesh for FEPE_SYN model runs. 

c = A*f,or A = c/f [m] (2-1) 
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where c is sound speed [m/s] , f is acoustic frequency [Hz] and 
X is wavelength [m] . The maximum frequency and minimum sound 

speed determine the smallest X. Therefore, the range/depth 

mesh was iteratively determined to be dr = A/2 = 0.206 m and 

dz = A/iO = 0.041 m, where fmax = 3600 Hz, c min = 1484 m/s 

for the data listed in Table 1. 

This range/depth mesh makes the computational time 
tremendously long, and a typical model run takes more than 100 
hours for the first of three phases of the computation. 
However, the objective of this research is to model energy 
spreading in the time domain, and measure TL and ESL 
quantitatively. Therefore, the long model run times are 
necessary and the development of a faster and more efficient 
program code is left for future research. The implementation 
of FEPE and FEPE_SYN on a massively parallel i860 array 
processor in a VMEbus is highly recommended for future ESL 
modeling computations at tactical sonar frequencies. 

B. NEW QUANTITATIVE MEASURE OF ESL 

Three measures to quantify ESL were introduced in the 
previous chapter (time spreading o, Jones ESL and MML) . 
Another measure of ESL is presented in this section and is 
based on the ability to model the time stretching of the 
signal accurately in shallow water. 

ESL is the result of time stretching of the transmitted 
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pulse due to multipath propagation in shallow water. The peak 
energy of the received pulse is reduced due to this time 
spreading. To measure ESL, consider a pulse at some range from 
the source which has undergone no spreading and which is 
symmetric in amplitude about its peak value in the time 
domain. This idealized “no spread” pulse at range, r, has the 
same time duration as the transmitted pulse at range r = l m. 
The quantitative measure of ESL is based on the pulse source 
level or the total energy at 1 m from the source: 



where E represents the total energy of the pulse, A represents 
the amplitude of the pulse at time t with a time increment dt 
and k is the number of time increments in the pulse, typically 
16,384 in this study. Erepuca is used here to compare its 
magnitude to the received pulses. 

In defining ESL quantitatively, an idealized, fictitious 
"no spread" pulse was determined by time compressing all of the 
pulse's time stretched energy into a pulse with no spreading 
that has the same pulse length and amplitude shape as the 
original transmitted Blackman pulse (Erepuca) • 

Assuming that the shape of both the replica and the no 
spread pulses are symmetric in amplitude, the ratio is 
defined: 




( 2 - 2 ) 
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(k) . 



(2-3) 



— ^Viospread (k) / ^replica 



Since energy must be conserved at the target within the time 
window, for both the idealized “no spread" pulse and the actual 
time spread pulse, require 



■•nospread 



= E 



Spread 



. Ea= 



spread 



(k) 



(2-4) 



Substituting Equations (2-2) and (2-3) into Equation (2-4) 
yields 



-•nospread 



= E 



^ nospread (k) — r * ^replic 



(2-5) 



Thus, the factor r is given by: 



r = (E, 



nospread 



/ E 



1/2 



replica 



( 2 - 6 ) 



Therefore, one can obtain the peak value, A^ospread'"^’' (k) ^ from the 
known values of A,-epiica'"®Mk) ^ E ^epiica and E spread as determined 
from Equation (2-4). Figure 1 shows the relationship of the 
peak values of the replica, no spread, and spread pulses. 

By defining these quantities as given by Equations (2-2) 
through (2-6) above, ESL and TL are automatically separated. 

TL is the energy lost between the source and target and 
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is calculated by FEPE_SYN in the form: 

TL = 10*logio (A.ospread'^'VA.eplica ' 

+ 10*logio (R) [dB] (2-7) 

where R is the range [m] from the source. ESL is the energy 
level that reaches the target that is spread in time beyond 
that of the original transmitted pulse length. 

In order to determine how much energy is lost due to time 
spreading, the squared maximum amplitude of the spread pulse 
is compared to the squared maximum amplitude of the no spread 
pulse, and their ratio is defined as ESL: 

ESL = 10*lOgi0 (Anospread""’' / A,p,ead""* )' 

= 20*log:o (A^oapread'""’' / ) [dB] (2-8) 

where Anospread”"^* is the absolute value of the pealc amplitude 
(plus or minus) of the no spread pulse. Aspread"'^’' is the 
absolute value of the peak amplitude (plus or minus) of the 
largest time- stretched pulse. 

A unique feature of the approach in this study is that 
the spread time series is modeled using very accurate, 
computer intensive models. Hence, no noise or artifacts 
contaminate this computation of ESL. Also, unlike previous 
studies, an assumed statistical distribution (e.g., Gaussian) 
for the time spread signal is not necessary; it is modeled 
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directly . 
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Table 1. Typical FEPE_SYN input data file 



Input Variables 


Comments 


-200 800.0 .06103515 16384 


(TMIN, TMAX, DT, NDT) 


3500.0 7.30 10.30555555 


(FREQ, ZS, ZR) 


100 15000 


(RSTART, RFINISH) 


15000 .206111111 50 


(RMAX, DR, NDR) 


110 .041222222 25 


(ZMAX, DZ, NDZ) 


0.0 64.1 


(ZSTART, ZFINISH) 


6000 


(ISIZE) 


1484.0 2 


(CO, NPADE) 


1 1.0 89.0 


(ISTRT, RMIN, THMAX) 

<Profile block (Hamilton 
geoacoustic model) follows 
as below> 


O.OOOE+00 64.1 
-1.000 -1.000 


(RD, D) 


★ * 

★ ★ 

-1.000 -1.000 


(Z, CW) 


* * 

★ ★ 

-1.000 -1.000 


(Z, CB) 


★ ★ 

★ ★ 

-1.000 -1.000 


(Z, RHOB) 


★ ★ 

★ ★ 


(Z, ATTN) 


-1.000 -1.000 


<can repeat more profile 
blocks to model a range 
dependent model > 
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Amplitude 





Relative Time (sec) 



Relative Time (sec) 



Figure 1. Relationship of the replica, no spread and spread 
pulses. Left: Amplitude of the pulses in the time domain. 
Right: Power (or acoustic energy) of the pulses in the time 
domain. 
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III. MODEL RESULTS 



A. BACKGROUND 

The transmitted pulse used in this study is a Blackman 
windowed pulse centered at 3500 Hz with a 200 Hz bandwidth, as 
shown in Figure 2. The Blackman pulse in the time domain is 
modeled with a pulse duration of 0.02 sec or 20 msec and is 
approximated by the FFT subroutine in EXT_TD. The received 
pulses are the Fourier synthesized time domain products of 200 
transmitted single frequency signals, each subject to an ocean 
transfer function (similar to TL) . Typical examples 
illustrating varying degrees of time stretching for a 
negative sound speed profile, overlying a sand bottom are 
shown in Figures 3 through 5. Figure 3 represents a situation 
with limited time stretching (ESL = l ~ 2 dB) , Figure 4 
moderate time stretching (ESL = 3 ~ 5 dB) and Figure 5 

extensive spreading (ESL =9-11 dB) . 

Four SSPs were used to determine the nature of ESL under 
widely varying, but common, shallow water SSPs (Figure 6) . 
Here the negative type SSP (negative SSP) corresponds to the 
observational data of Area Foxtrot [Scanlon, 1995] . The 
isothermal type (isothermal SSP) represents well mixed 
conditions often found in winter in shallow water regions. The 
mixed layer type SSP (ML- type SSP) is representative of summer 
conditions where the water column has been warmed by 
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insolation but its upper portion has been mixed by wind 
forcing. The deep sea sound channel type SSP (SC- type SSP) is 
introduced to represent a typical deep water situation where 
little time stretching is anticipated at short ranges for a 
shallow source. Because of the extensive computer time 
required to perform the analysis for a typical 4000 m deep 
water column, the lower positive gradient of this profile has 
been increased to mimic deep water refraction patterns. 

In order to examine the impact of a slow speed sediment 
bottom on ESL, the isothermal profile was arbitrarily 
increased by 70 m/sec ( (b) in the upper right panel of Figure 
6) to artificially create a slow speed sediment interface 
(slow speed SSP) . Hamilton geoacoustic models for both sand 
and silt/clay bottoms were used as inputs to FEPE (Tables 2 
and 3) . These geoacoustic models correspond closely with the 
measured geological data in Area Foxtrot [Scanlon, 1995] . 

Because it is iirpractical to make three dimensional color 
plots for a pulse centered at 3.5 kHz with a 200 Hz bandwidth, 
two dimensional color plots at a single frequency, calculated 
by FEPE, are shown to illustrate propagation characteristics 
in shallow water. The horizontal and vertical resolution are 
10.31 m and 0.21 m, respectively, for all the FEPE TL color 
plots. Figures 7 and 8 depict the TL for a negative SSP 
overlying a sandy bottom and a clay/silt bottom, respectively. 
A comparison of these two figures illustrates the impact of 
sediment type on TL. Because of the high attenuation 
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associated with a silt/clay bottom, the TL is significantly 
greater than for a sand bottom. Figure 9 is a plot of the TL 
for an isothermal (weakly positive) SSP overlying a sand 
bottom, and Figure 10 is a plot of TL for a mixed layer type 
(ML- type) SSP, sand bottom showing the favorable propagation 
conditions within the ML. 

Figure 11 is an enlarged view of the TL shown in Figure 
8 for the negative SSP silt/clay bottom case but for 
frequencies of 3400 Hz and 3600 Hz, the minimum and maximum 
frequencies within the bandwidth. The dispersion of acoustic 
energy over the 2 00 Hz bandwidth is readily noted. A 
difference of up to a 100 m is observed where the bottom 
bounce rays reflect from the bottom near 5000 m range. When 
modeling broadband pulses, TL computed by FEPE can only be 
considered a gross approximation to the propagation of finite 
bandwidth pulses because of this dispersion effect, but it can 
illustrate general TL properties for various environmental 
effects (e.g., SSP, sediment types, etc.) at the center 
freq[uency of the transmitted pulse. However, it is important 
to realize that FEPE TL estimates are for a single frequency 
only, with no temporal information about the pulse 
distribution, so that it can not be used to estimate ESL from 
a 200 Hz band sonar. 

The ESL plots in the following sections were generated on 
workstation (UNIX) computers using MATLAB programs based on 
the definitions in the Chapter II. These MATLAB programs. 
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developed by the author, use the Cray/Onyx computer outputs of 
the time- stretched signal and calculate and display ESL based 
on Equation (2-8). 

B. RANGE DEPENDENCE OF ESL 

Figure 12 shows ESL as a function of range between 0.1 km 
and 7.6 km for a negative SSP over a sand bottom. This plot 
has a range resolution of 10.3 m and ESL is plotted for 
targets positioned at 10.3 m depth intervals, i.e., from 10.3 
m to 61.8 m. The range dependence of ESL for the transmitted 
pulse is not a linear function of range as might intuitively 
be expected, i.e., greater ranges would imply more time 
stretching of the pulse. However, as Figure 12 shows, ESL 
increases significantly for ranges from 0 m to 1600 m, but 
remains relatively constant at longer distances as the 
transmitted pulse begins to saturate in its interaction with 
the shallow water boundaries . As the energy propagates farther 
in range, the higher order normal modes begin to separate from 
the lower normal modes due to their different group speeds. 
However, the higher order mode energy (steeper ray equivalent 
angles) is continually attenuated at a relatively faster rate 
due to bottom boundary interactions. Beyond 1600 m most of the 
energy is carried by the lower order modes with little 
variation in group speed. Hence, the ESL remains fairly 
constant over this distance. This plot also demonstrates how 
ESL varies with depth of the target. ESL remains relatively 
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high at all depths (up to a maximum of 11 dB) over this entire 
range band, but is higher in the mid-depths (20.6 m to 41.2 ra) 
due to the interference pattern of the different normal modes. 
Table 4 presents the mean and standard deviation of ESL for 
each depth over the range band 2000 m to 14500 m calculated at 
10.31 m intervals. Over this range band the mean ESL varies 
from approximately 6.5 to 7.5 dB with a standard deviation of 
about 1.5 dB . 

C. DEPTH DEPENDENCE OF ESL 

A plot of ESL versus depth is shown in Figure 13, which 
is based on a negative SSP over a sand bottom. These plots, 
based on a 1.031 m vertical resolution, pictorially describe 
how ESL varies throughout the water column at the ranges 
indicated above the individual panels (approximately 1000 m 
intervals) . A weak depth dependence of ESL for depths above 15 
m is noted due to the location of the source at 7.3 m. ESL is 
lower by 5 dB or more in the upper 15 m at selected ranges 
where the downward refracted energy from the shallow source 
returns to the near surface. In general, ESL is relatively 
independent of depth as noted in the previous section. Both 
high and low frequency fluctuations of 1 to 3 dB amplitude 
occur throughout the water column. For this source/SSP 
configuration strong high frequency fluctuations (1 to 6 dB) 
are observed between 20 m and 40 m. 

The near constancy of ESL with depth is a result of the 
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coherent suinination of approximately 290 to 307 modes 
propagating at different angles and at different group speeds 
at each frequency in the 200 Hz-wide band. Although higher 
order modes are attenuated faster, the interference pattern of 
the remaining modes results in smoothing the energy 
distribution throughout the water column. The fluctuations 
represent specific depth/range positions where constructive or 
destructive interference occurs. 

This relative lack of depth dependence for ESL appears to 
be fairly universal as investigations were done for a variety 
of SSP shapes and bottom sediment conditions (discussed later) 
and all cases appeared to support this finding. 



D. GEOACOUSTIC DEPENDENCE OF ESL 

Significant changes occur in TL when the geoacoustic 
parameters and sediment properties exhibit high spatial 
variability in shallow water environments [Scanlon, 1995] . For 
a location where the sediment compressional sound speed 
decreases due to the bottom transitioning from sand to 
silt/clay, one expects the received pulses to become fewer in 
number and show a significantly increased TL. Figure 14 shows 
ESL vs depth plots, for the isothermal SSP overlying a 
silt/clay bottom sediment. However, to model a “slow speed 
bottom", as described previously, the SSP in the water column 
has artificially been increased by 70 m/sec to create a low 
speed sedimentary layer. Because of the highly absorptive 
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nature of this slow speed bottom, only a few low order modes 
of similar group speed dominate the water column. Hence the 
number of time stretched pulses is expected to decrease. Only 
one pulse was observed for this case compared to 10 or more 
for the same SSP overlying a sand bottom. As seen in Figure 
14, there is virtually no ESL observed at any range or depth. 

As the bottom becomes more reflective, the amount of ESL 
is expected to increase. This is borne out by Figures 15 and 
16 which show ESL vs depth plots, for an isothermal SSP 
overlying a silt/clay and a sand bottom, respectively. 
Comparing Figure 14 (slow speed bottom) , Figure 15 (clay) and 
Figure 16 (sand) , mean ESLs averaged over the depth of the 
water column are seen to increase from 0.1 dB, 2.7 dB, 6.5 dB, 
respectively. Although these relative changes in ESL were 
based on the slightly positive isothermal SSP, the dependence 
of the water column SSP is relatively weak as similar ESL 
values were obtained when a negative SSP was substituted. This 
dependence of ESL on the water column SSP is examined in 
further detail below. 

E. SSP DEPENDENCE OF ESL 

The influence of the shape of the SSP of the water column 
on ESL can be examined by varying the SSP but keeping the 
geoacoustic properties of the bottom constant. Three SSPs were 
considered: the ML-type SSP (Figure 17), the negative SSP 
(Figure 13) and the isothermal SSP (Figure 16), all overlying 
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a sand bottom. The ML- type SSP is essentially a combination of 
the isothermal and the negative SSP. As can be seen in Figure 
17, ESL remains moderately low (~ 6.5 dB) in the mixed layer, 
but becomes larger (~ 9 dB) below the mixed layer (ML) with 
large scale fluctuations (2 ~ 3 dB) frequently occuring. In 
fact, these below layer ESL values were the highest 
encountered of all SSP/bottom type configurations examined. 
For the negative SSP (Figure 13) , acoustic energy is forced 
towards the bottom, so that ESL becomes large (average of 7.0 
dB for each depth column beyond 5152 m) . For the isothermal 
SSP (Figure 16) , acoustic energy is forced upward weakly, 
minimizing the amount of energy which penetrates into the 
bottom. For the reflective sand bottom a relatively moderate 
ESL is observed (~ 6.5 dB) for ranges beyond 5152 m. 

These three examples demonstrate that the shape of the 
water column SSP exerts a relatively insignificant influence 
on ESL. For the reflective bottom type considered, all three 
yielded nearly similar values (7 ~ 8 dB) . Negative SSPs cause 
ESL values to be slightly greater (- 1 dB) than upward 

refracting SSPs. This variability is small compared to the 3 
~ 5 dB difference noted between an absorptive (silty) and 
reflective (sandy) bottom. 

Because the SSP shape has only a limited effect on ESL, 
this suggests there will be only a limited seasonal 
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variability associated with ESL for any given shallow water 
location. 
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Table 2. A Hamilton geoacoustic model for a sand bottom in 
Area Foxtrot 



depth 



0 


OOOOOOE+00 


1519 


000 


6 


250000 


1518 


000 


12 


50000 


1516 


. 000 


18 


75000 


1512 


. 000 


25 


OOOOO 


1509 


000 


31 


25000 


1506 


, 000 


37 


50000 


1500 


000 


43 


75000 


1496 


000 


50 


OOOOO 


1493 


000 


56 


25000 


1490 


000 


62 


50000 


1488 


000 


64 


OOOOO 


1484 


000 


64 


. 10000 


1484 


000 


1 


000000 


- 1 . 000000 


64 


10000 


1817 


0 0 0'' 


65 


ooooo 


1934 


400 


66 


ooooo 


1934 


, 850 


67 


ooooo 


1935 


290 


68 


ooooo 


1935 


720 


69 


ooooo 


1936 


140 


70 


ooooo 


1936 


560 


71 


ooooo 


1936 


970 


”7 2 


ooooo 


1937 


380 


''3 


ooooo 


1937 


780 


74 


ooooo 


1938 


180 


75 


ooooo 


1938 


570 


76 


ooooo 


1938 


950 


7 7 


ooooo 


1939 


330 


78 


ooooo 


1939 . 


710 


79 


ooooo 


1940 


080 


80 


ooooo 


1940 . 


440 


81 


ooooo 


1940 . 


810 


82 


ooooo 


1941 


000 


83 


ooooo 


1941 


520 


64 


ooooo 


1941 


870 


85 


ooooo 


1942 


210 


86 


ooooo 


1942 


550 


87 


ooooo 


1942 


890 


88 


ooooo 


1943 


220 


89 


ooooo 


1 943 . 


550 


90 


ooooo 


1943 


880 


91 


ooooo 


1944 


200 


92 


ooooo 


1944 


520 


93 


ooooo 


1944 


830 


94 


ooooo 


1945 


140 


94 


10000 


2000 


000 


lie 


1 0 0 0 0 


2000 


000 


1 000000 


1 000000 


6 4 


10000 


2 600000“ 


65 


ooooo 


2 600000 


b i 


ooooo 


2 600000 


6 7 


ooooo 


2 610000 


6 6 


ooooo 


2 610000 


69 


0000(- 


2 620000 


70 


0000c 


2 6 20 0 0 0’ 


i 


0000 (■ 


: b 1 0 0 0 (■ 


*7 2 


ooooo 


: 6 3 0 C' 0 (• 



73 00000 

74 00000 

75 00000 

76 00000 

77 00000 

78 00000 

79 00000 

80 00000 
81 00000 
82 00000 

83 00000 

84 00000 

85 00000 

86 00000 

87 00000 

88 00000 

89 00000 

90 00000 

91 00000 

92 00000 

93 00000 

94 00000 
94 10000 
110 0000 
1 000000 

64 10000 

65 00000 

66 00000 

67 00000 

68 00000 
69 00000 

compressional 7o ooooo 
wave speed 7i ooooo 
(m/s) ooooo 

73 ooooo 

74 ooooo 

75 ooooo 
"6 ooooo 

77 ooooo 

78 ooooo 

79 ooooo 

80 ooooo 

81 ooooo 

82 OOOOO 

83 OOOOO 

84 OOOOO 

85 OOOOO 

86 OOOOO 

87 00000 

88 OOOOO 
69 OOOOO 
90 OOOOO 
<^1 OOOOO 

92 00000 

93 OOOOO 

>^4 OOOOO 
“-I 10000 

; j 0 0 0 0 

. 000000 



water column 
sound speed 
(m/s) 



2 640000 
2 640000 
2 650000 
2 650000 
2 650000 
2 660000 
2 660000 
2 670000 
2 670000 
2 680000 
2 680000 
2 690000 
2 690000 
2 720000 
2 740000 
2 760000 
2 780000 



sediment 

density 

(g/cm^) 



2 800000 
2 810000 
2 820000 
2 830000 
2840000 



2 . 850000 
2 850000_ 
- 1 000000 
0 703474 0~ 
0 8056971 
0 8307893 
0 . 8309792 
0 8311670 
0 8313505 
0 8315299 
0 8317071 
0 8318844 
0 8320595 
0 8322304 
0 8323971 
0 8325637 
0 8327284 
0 8328909 
0 8330534 
0 8332115 
0 8333656 
0 8335091 
0.8336610 
08338214 
0 8339733 
0 8341210 
0 8342644 
0 8344080 
0 8345492 
0 8346906 
0 8348320 
0 8349692 
0 8351042 
0 8352392 
0 8396534 | 
20 OOOOO J 



sediment 
attenuation 
(dB/A ) 



1 000000 
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Table 3. A Hamilton geoacoustic model for a silt/clay 
bottom in Area Foxtrot 



depth 

0 . OOOOOOE+OO 


1519.000 




73 00000 

74 00000 

75 00000 


1 . 540000 
1 540000 
1 . 550000 


6 . 250000 


1518 . 000 




76 00000 


1 . 550000 


12 . 50000 


1516 . 000 




77 00000 


1 550000 


18 . 75000 


1512.000 




78 00000 


1 . 560000 


25 . 00000 


1509 000 


water column 


79 . 00000 


1 . 560000 


31 . 25000 


1506 . 000 


80 . 00000 


1 . 570000 


37 . 50000 


1500.000 


sound speed 


81 . 00000 


1 . 570000 


4 3 .75000 


1496 . 000 


(m/s) 


82 .00000 


1 . 580000 


50 . 00000 


1493 . 000 


83 00000 


1 580000 


56 . 25000 


1490 . 000 




84 . 00000 


1 . 590000 


62 . 50000 


1488 . 000 




85 . 00000 


1 . 590000 


64 .00000 


1484 . 000 




86 . 00000 


1 620000 


64 . 10000 


1484 . 000 _ 




87 00000 


1 640000 


- 1 . 000000 


- 1 . 000000 




88 .00000 


1 660000 


64 . 10000 


1550 . 000 “ 




89 .00000 


1 .680000 


65 00000 


1550 . 500 




90 . 00000 


1 .700000 


66 .00000 


1555 . 000 




91 . 00000 


1 .710000 


67 . 00000 


1560 . 000 




92 . 00000 


1 .720000 


68.00000 


1565.000 




93 . 00000 


1 730000 


69 . 00000 


1570 . 000 




94 . 00000 


1 .740000 


70 00000 


1580 100 




94 10000 


1 .750000 


71 . 00000 


1584 .000 




110 . 0000 


1 . 750000 _ 


72 .00000 


1587 . 300 




-1 000000 


- 1 . 000000 


73 . 00000 


1590 300 




64 . 10000 


0 . 400000 “ 


74 00000 


1593 .000 




65 00000 


0 410000 


75 . 00000 


1595. 400 




66 . 00000 


0 420000 


76 .00000 


1597 . 600 


compressional 


67 . 00000 


0 .430000 


77 . 00000 


1599 . 700 


wave speed 
(m/s) 


68 . 00000 


0 .440000 


78 . 00000 


1601 . 600 


69 . 00000 


0 450000 


79.00000 


1603 . 300 


70 . 00000 


0 .460000 


80 00000 


1605 . 000 




71 . 00000 


0 470000 


81 00000 


1606 . 500 




72 00000 


0 480000 


82 .00000 


1608 . 000 




73 .00000 


0 490000 


83 00000 


1609.400 




74 . 00000 


0 500000 


84 00000 


1610 .700 




75 00000 


0 510000 


85 . 00000 


1611 . 900 




76 00000 


0 520000 


86 . 00000 


1613 100 




77 00000 


0 . 530000 


87 .00000 


1614 300 




78 .00000 


0 . 540000 


88 .00000 


1615 400 




79 .00000 


0 . 550000 


89 . 00000 


1616 . 400 




80 . 00000 


0 . 560000 


90 . 00000 


1617 . 400 




81 . 00000 


0 . 570000 


91 . 00000 


1618 . 400 




82 . 00000 


0 . 580000 


92 . 00000 


1619 . 400 




83 00000 


0 . 590000 


93 . 00000 


1620 300 




84 . 00000 


0 . 600000 


94 00000 


1621 100 




85 . 00000 


0 .610000 


94 10000 


1621 . 200 




86 00000 


0 .620000 


110 0000 


1632 200 ^ 




87.00000 


0.630000 


1 000000 


1 000000 “ 




88 00000 


0 640000 


64 10000 


1 500000 “ 




89 00000 


0 650000 


65 .00000 


1 500000 




90 .00000 


0 . 660000 


66 00000 


1 500000 




91 00000 


0 670000 


67 00000 


1 510000 




92 00000 


0 680000 


68 00000 


1 510000 




93 00000 


0 690000 


69 00000 


1 520000 




94 00000 


0.700000 


70 00000 


1 520000 




94 10000 


0 710000 


71 00000 


1 530000 




110 0000 


20 . 00000 _ 


72 00000 


1 530000 




• 1 000000 


1 000000 ” 



sediment 

density 

(g/cm^) 



sediment 
attenuation 
(dB/X ) 
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Table 4. The mean and standard deviation of ESL 



Depth (m) 


Mean ESL(dB) 


Std ESL(dB) 


10.3 


5 . 0364 


1.712 


20.6 


6.9412 


1.474 


30.9 


7.2397 


1.466 


41.2 


7.6396 


1.429 


51.5 


6.7190 


1.210 


61.8 


5.9725 


1.171 
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Normalized Power Norrrialized Amplitude 
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Figure 2. A Blackman windowed pulse, used as the input 
transmitted pulse in this research, is shown in the 
frequency domain (left) and in the time domain (right) . The 
upper panels are the amplitudes normalized by the peaks; 
the lower panels are the power (or acoustic energy) 
normalized by the peaks. 
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Normalized Amplitude 





Figure 3 . An example of time stretching illustrating low 

ESL (~ 2 dB) for a negative SSP overlying a sand bottom; 

source depth is 7.3 m. 
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Normalized Amplitude 




Figure 4. Same as Figure 3 except illustrating moderate ESL 
{- 4 dB) . 
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Figure 5. Same as Figure 3 except illustrating large ESL (~ 
10 dB) . 
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Depth (m) Depth (m) 




Sound Speed (m/sec) 




Sound Speed (m/sec) 




Sound Speed (m/sec) 




Sound Speed (m/sec) 



Figure 6. Sound speed profiles (SSP) used in this research: 
"negative SSP” (upper left) , "isothermal SSP” (upper right) , 
“mixed layer (ML) - type SSP” (lower left), "sound channel (SC) - 
type SSP" (lower right) . 
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Figure 7. Three dimensional TL plots by FEPE for 3.5 kHz 
(single frequency) , for a negative SSP overlying a sand 
bottom; water depth is 64 m, source depth is 7.3 m. 
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Figure 8. Same as Figure 7 except for a silt/clay bottom. 
The color bar scale was changed due to large TL for this 
case. 
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Figure 9 , Same as Figure 7 except an isothermal SSP is 
used. 
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Figure 10. Same as Figure 7 except a mixed layer type SSP 
is used. 



43 








i 



i* 



Depth (m) Depth (m) 



Frequency = 3400 Hz 





4000 



4500 5000 

Range (m) 

Frequency = 3600 Hz 



5500 



6000 TL (dB) 



Or 

10 - 

20 - 

30- 

40- 

50- 

60- 



4000 4500 5000 5500 6000 TL (dB) 

Range (m) 



Figure 11. Three dimensional TL plots by FEPE for 3.4 kHz 
(upper) and 3.6 kHz (lower) (single frequency), for a 
negative SSP overlying a silt/clay bottom; water depth is 
64 m, source depth is 7.3 m. 
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Figure 12. ESL versus range plots for a negative SSP 
overlying a sand bottom. Source depth is 7.3 m. Target 
depths are 10.3, 20.6, 30.9, 41.2, 51.5, 61.8 m, 
respectively . 
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Figure 13 . ESL versus depth plots for a negative SSP 
overlying a sand bottom. Source depth is 7.3 m. Ranges 
shown on the top of each panel . 
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Figure 14. ESL versus depth plots for an isothermal SSP 
overlying a slow speed bottom. Source depth is 7.3 m. 
Ranges are shown on the top of each panel . 
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Figure 15. Same as Figure 14 except for a silt/clay bottom. 
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Figure 16. Same as Figure 14 except for a sand bottom. 
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Figure 17. ESL versus depth plots for a ML- type SSP 
overlying a sand bottom. Source depth is 7.3 m. Ranges 
shown on the top of each panel . 
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IV. OTHER ESL MEASURES 



As discussed previously, several techniques have been 
introduced by various researchers to measure the amount of 
time stretching caused by multipath propagation. All suffer 
some degree of degradation due to assumptions inherent with 
each technique. In this chapter a quantitative comparison of 
three commonly used techniques is made with the ESL technique 
developed in this study. 

A. TIKE SPREADING STANDARD DEVIATION 

The standard deviation of the time spreading, o, (Equation 
(1-2)) and the average arrival time (Equation (1-3)) were 
calculated by a MATLAB program developed by the author for all 
data processed. Figures 18 and 19 show the normalized (by the 
peak value) ESL and o as a function of depth and range, 
respectively. These two very different parameters compare 
favorably in a qualitative sense. However, there are depths 
and ranges where the correlation between ESL and o is not 
good, and may be due to the fact that the distribution of the 
pulse time spreading is definitely not Gaussian (see Figures 
3 through 5) . Hence, the standard deviation may not be as 
meaningful as it would be if the distribution of the energy 
within the pulse were Gaussian. 
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B. JONES' DEFINITION OF ESL 

ESL based upon the definition of Jones (1990) (Equation 
(1-4)) was also computed by a MATLAB program developed by the 
author. Jones' theoretical approach is similar to the 
definition of ESL given in this study for a situation of a 
signal uncontaminated by noise. Both Equations (1-4) and (2-8) 
will yield identical results if the size of the resolution 
cell (the temporal resolution of the target system) is 
selected accurately. Unfortunately, it is difficult to 
determine the optimum duration of the resolution cell without 
first modeling ESL accurately. Jones suggested (1990) using a 
resolution cell width ec[uivalent to the reciprocal of the 
pulse bandwidth, for this study 1/200 = 0.005 sec. 

When the pealc pulse amplitudes of the time stretched 
pulses from the two techniques are compared, it is found that 
the Jones technique overestimates the magnitude of ESL by 2.4 
dB. In order to bring the Jones technique into agreement with 
that of this study, the 2.4 dB reduction in peak amplitude is 
equivalent to increasing the resolution cell size to 0.0084 
sec. The 0.0084 sec resolution cell size is not a general 
solution but is appropriate for short duration active pulses 
near 3.5 kHz. Accurate modeling of the time- stretched pulse is 
required to overcome the potential error in ESL if the 
resolution cell size is determined by a best guess. 
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C. MML 

A normalized cross correlation of Equation (1-6) was also 
calculated by a MATLAB program developed by the author. 
Because of the simple temporal pulse shape, the coherence 
values were very high near the peak of the propagated pulse 
(Table 5) with a maximum value of 0.9999 obtained for the 
isothermal SSP, silt/clay case. This result supports the 
validity of the assumption that the shape of the pulse with no 
spreading and no reverberation is symmetric to the transmitted 
pulse described in Section B of Chapter II. However, the 
received pulse shape is not perfectly symmetric to the 
transmitted one. It is a little thinner than expected because 
the normalized cross correlation never reaches unity in this 
research, and several ESL values of around -0.0029 dB were 
observed for the isothermal SSP, silt/clay case. In the real 
ocean environment the MML will be severely degraded by the 
presence of reverberation (and ambient noise) , and the high 
MML values obtained here for the model/signal-only case (i.e., 
no noise contamination) are expected to be high. A correlation 
beamformer with an advanced signal processor, such as Inverse 
Beamforming (IBF) , offers great promise in regaining/ 
overcoming sonar system performance degradation due to ESL in 
shallow water. 

For the tactical active sonar, matched filter or 
correlation processing is the heart of the detection system. 
When the signal gives low coherence (equivalent to the 
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normalized cross correlation given by Equation (1-6)) below a 
previously selected threshold, it is recognized as a false 
signal. Jensen and Sabbadini (1993) demonstrated the impact of 
the signal coherence by defining the MML for a LFM signal with 
bottom interaction in deep water. One expects this may have 
significance when sophisticated pulses such as LFM signals are 
used in shallow water. Coherenced- based signal processing/ 
beamforming methods used with pea)c picJcers for post processors 
show great promise in regaining ESL due to time 
stretching. [Nuttall and Wilson, 1991; Wilson, 1995; Fabre and 
Wilson, 1995] 
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Table 5. The coherence (normalized cross correlation) for an 
isothermal SSP overlying a silt/clay bottom 



Range (m) 


Mean 

coherence 


std 

coherence 


Max 
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Min 

coherence 
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. 7715 
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Figure 18. Comparison of ESL to time spreading standard 
deviation for a negative SSP overlying a sand bottom. 
Source depth is 7.3 m, target depth is 10.3 m. The solid 
line represents ESL, the dots represent the standard 
deviation . 
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Figure 19. Comparison of ESL to time spreading standard 
deviation for a negative SSP overlying a sand bottom. 
Source depth is 7.3 m, range is 2504 m. The solid line 
represents ESL, the dots represent the standard deviation. 



57 




58 



V. IMPACT OP ESL ON TACTICAL ACTIVE SONAR PERFORMANCE IN 
SHALLOW WATER 

A. IMPACT OF ESL ON TOTAL LOSS 

According to the ESL analysis in the previous chapter, it 
is the bottom boundary geoacoustic properties that dominate 
the behavior of both TL and ESL. Thus, a tactically important 
SSP is one exhibiting a negative profile because this SSP 
refracts acoustic energy downward towards the bottom boundary. 
It was the negative SSP overlying sand (reflective) bottom 
which exhibited the highest ESL. In contrast, when this SSP 
overlies a silt/clay (absorptive) bottom, the ESL is 
significantly reduced (4 dB or more) , but as expected, the TL 
is also very large. 

Figure 20 shows an example of the received pulses for the 
negative SSP overlying a silt/clay bottom. Compared to the 
stretching observed in Figures 3 through 5, one sees for this 
case that the time stretching is minimal; the pulses 
associated with the higher modes which lead to the time 
dispersion are absorbed by the bottom. This is borne out in 
the ESL vs range plots where the sand (reflective) bottom 
(Figure 21) , shows ESL to be ~ 3 to 4 dB higher than for the 
silt/clay (absorptive) bottom (Figure 22) . 

In spite of the low ESL for a silt/clay bottom, the 
overall detection range is degraded because TL is extremely 
large. Thus, a reflective sand bottom probably offers the best 
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detection opportunity because of the reduced TL. If signal 
processing methods can be developed to regain the loss due to 
energy spreading, then sandy, reflective bottoms can be 
anticipated to offer even greater detection opportunities. 

The impact of the bottom sediment on sonar performance is 
further illustrated by comparing the total transmission loss 
(ESL + TL) and TL only for propagation over sandy and silty 
bottoms (Figures 23 and 24) . The difference between the curves 
on each plot is the ESL. Note that as stated previously, the 
two curves are almost coincident for the silt/clay bottom 
(i.e., ESL is small) but widely separated for the sand bottom 
(ESL large) . However, better sonar performance is not realized 
in areas covered by silt/clay (absorptive) sediments because 
the potential gain achieved by low ESL is offset by a much 
larger TL. 



B. ANALYSIS OF ESL IN DEEP WATER 

Although ESL is significant in shallow water environments 
due to multipath effects, ESL isn't as large in deep water 
because a single, or nearly similar multipaths, dominates the 
time arrival structure. For passive sonars at very long ranges 
one can measure significant time dispersion even in deep 
water, but tactical active sonar detection ranges are usually 
very short. Ray theory is usually a good approximation and the 
time stretching is usually minimal in deep water. Figure 25 
shows the TL for a single frequency 3500 Hz signal in shallow 
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water where a deep sound channel (SC) sound speed profile with 
axis at 25 m has been simulated to illustrate this ESL deep 
water dependence. 

Although this is not a realistic situation, Figure 25 
shows the impact of the SC- type SSP on TL. ESL vs depth plots 
for the received pulses are shown in Figure 26. Comparing the 
impact of ESL to the negative SSP case (Figure 13) , ESL is far 
less significant for the deep water, SC propagation. Figure 27 
shows TL + ESL vs depth plots for this case. It is evident 
that at depths where TL is small, ESL is also small. Therefore 
ESL is not significant in the deep water, short range SC 
propagation environment. 



C. ANALYSIS OF DEEPER SOURCE DEPTH 

So far the analysis has been performed only for the 
source depth (SD) of 7.3 m, chosen to represent the typical 
depth for a bow mounted sonar. However, we must consider 
variable depth, tactical sonars which can be deployed from an 
ASW helicopter. Although this sonar can change its operational 
depth, the optimum depth for performance can be found by 
computing a nximber of FEPE runs at various source depths for 
a single frequency (3500 Hz) . We selected a source depth of 
51.5 m as optimum based on analysis. Figures 28 and 29 depict 
the TL for a negative SSP overlying a sand and silt/clay 
bottom, respectively. Figure 30 and 31 show the ESL vs depth 
plots for the above two cases. Comparing Figure 13 (shallow 
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source) to Figure 30 for the sand bottom, ESL is seen to 
decrease as the source depth increases. Examining both 
received pulses in the time domain, the observed number of 
time spread pulses significantly decreased from the order of 
20 for the shallow source to the order of 5 for the deep 
source. This indicates that bottom interaction was decreased 
for a near bottom source, resulting in fewer pulses and less 
time spreading. The same result is obtained for the silt/clay 
bottom as seen in the comparison of Figures 22 and 31. In this 
case, even fewer pulses were observed for the deep source 
depth . 

In Figures 28 and 29, the direct path energy from a deep 
source propagates to great ranges (up to 8000 m range) . 
Accordingly, for the negative SSP with the deeper source 
depth, propagating energy for depths near the bottom is not 
attenuated as readily as for the shallow source. Propagating 
modes exist in the deeper depth region and, combined with the 
reduced ESL impact for the deeper source, implies that 
improved performance may be anticipated for the variable depth 
sonar when lowered to depths well below the surface. 

D. COMPARISON OP BROADBAND PULSE TO SINGLE FREQUENCY 

PROPAGATION 

Although FEPE is being adopted as the new Oceanographic 
and Atmospheric Master Library (OAML) standard acoustic 
propagation model by the U.S. NAVY, FEPE predicts the TL only 
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for a single frequency, which contains no temporal spreading 
information. In reality, pulses with finite bandwidths are 
used in all tactical active sonars. To illustrate the 
difference between a single and multiple- frequency TL, TL 
plots for both a single frequency (3.5 kHz) and a 200 Hz band 
pulse centered at 3.5 kHz are shown in Figure 32 for a 10.31 
m target depth for the negative SSP, sand bottom case. The 
thin line shows the single frequency (3.5 kHz) FEPE TL 
estimate, and the thick line shows the FEPE_SYN TL estimate 
for the broadband pulse. The 200 Hz pulse shows less TL 
fluctuations with range than the single frequency TL. This is 
the result of convolution (or Fourier synthesizing) of 200 
frequency bins of energy which behaves like incoherent 
summation. This is tactically important because for the 
shallow water wave guide, TL will not fluctuate severely, but 
will increase rather smoothly as will the probability of 
detection curve. To illustrate this effect on ESL Figure 33 
shows the total loss (200 Hz band) and single frequency TL 
(3.5 kHz) vs range for the same inpit parameters. ESL degrades 
the total loss by 1 to 8 dB. This means that a single 
frequency FEPE TL estimate can underestimate the total loss 
from 1 to 15 dB in shallow water. 

Consequently, for tactical active sonars, the performance 
should be predicted using the transmitted pulse TL, not the 
single frequency TL. 



63 



E. EFFECT ON SMALL SCALE CHANGES IN TARGET DISPLACEMENT 

The total loss vs range or depth plots (e.g., lower 
panels of Figure 23 and 27) indicate that small range and 
depth scale changes can possibly be significant. Because of 
computational efficiency, the resolutions of all the previous 
plots were selected with a range increment of 10.31 m and a 
depth increment of 1.03 m. Figure 34 shows ESL calculated for 
range increments, 10.31 m and 0.412 m for the negative SSP, 
sand bottom case. This figure shows that the coarse horizontal 
resolution replicates the features of the transmitted pulse 
satisfactorily. Similarly, a coarse (1.031 m) and a fine 
(0.206 m) depth increment intercompared . Figure 35 shows that 
the coarse vertical increment replicates the transmitted pulse 
satisfactorily but exhibits finer- scale fluctuations that are 
not reproduced by the coarse depth resolution. 

These figures also demonstrate that the total loss is not 
as sensitive to horizontal array displacements as for vertical 
displacements. Thus, the sonar array should be vertically 
stabilized physically or electronically. 



64 



Normalized Amplitude 




Figure 20. An example of the received pulse shape for a 
negative SSP overlying a silt/clay bottom, source depth is 
7.3 m. ESL is 2.0 dB . 
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Figure 21. ESL versus range plots for a negative SSP 
overlying a sand (reflective) bottom. Source depth is 7.3 
m. Target depths are 10.3, 20.6, 30.9, 41.2, 51.5, 61.8 m, 
respectively . 



66 



ESL (dB) 




Figure 22. Same as Figure 21 except for a silt/clay 
(absorptive) bottom. 
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Figure 23 . Total loss and TL versus range plots for a 
negative SSP overlying a sand (reflective) bottom. Source 
depth is 7.3 m. Target depths are 10.3, 20.6, 30.9, 41.2, 
51.5, 61.8 m, respectively. The thick line represents TL, 
the thin line represents total loss (TL + ESL) . 
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Figure 24. Total loss, TL versus range plots for a negative 
SSP overlying silt/clay (absorptive) bottom. Source depth 
is 7.3 m. Target depths are 10.3, 20.6, 30.9, 41.2, 51.5, 
61.8 m from the top panel, respectively. The thick line 
represents TL, the thin line represents total loss (TL + 
ESL) . 
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Figure 25. Three dimensional TL plots by FEPE for 3.5 kHz 
(single frequency) for a SC- type SSP overlying a sand 
bottom; water depth is 64 m, source depth is 7.3 m. 
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Figure 26. ESL versus depth plots for a SC- type SSP 
overlying a sand bottom. Source depth is 7.3 m. Ranges 
shown on the top of each panel . 
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Figure 27. Total loss and TL versus depth plots for a 
type SSP overlying a sand bottom. Source depth is 7.3 
Ranges are shown on the top of each panel . 
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Figure 28. Three dimensional TL plots by FEPE for 3.5 kHz 
(single frequency) for a negative SSP overlying a sand 
bottom; water depth is 64 m, source depth is 51.5 m. 
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Figure 29. Three dimensional TL plots by FEPE for 3.5 kHz 
(single frequency), for a negative SSP overlying silt/clay 
bottom, water depth is 64 m, source depth is 51.5 m. The 
color bar scale was changed due to large TL for this case. 
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Figure 30. ESL versus depth plots (upper panels) and total 
loss and TL versus depth plots (lower panels) for a 
negative SSP overlying a sand bottom. Source depth is 51.5 
m. Ranges are shown on the top of each panel . 
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Figure 31. Same as Figure 



30 except for a silt/clay bottom. 
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Figure 32. Comparison of TL from a single frequency (thin 
line) and TL from a 200 Hz band pulse (thick line) for a 
negative SSP overlying a sand bottom. Source depth is 7.3 
m, target depth is 10.3 m. 
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Figure 33. Comparison of TL from a single frequency (thin 
line) and total loss (TL + ESL) from a 200 Hz band pulse 
(thick line) for a negative SSP overlying a sand bottom. 
Source depth is 7.3 m, target depth is 10.3 m. TL from a 
200 Hz band pulse (dashed line) is also shown for 
reference . 
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Figure 34. Comparison of ESL at 10.31 m range inc 
(stars) and ESL at 0.412 m range increment (solid 1 in ) 
a negative SSP overlying a sand bottom. Source depth 
m, target depth is 10.3 m. 
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Figure 35. Comparison of ESL at 1.031 m depth increment 
(stars) and ESL at 0.206 m depth increment (solid line) for 
a negative SSP overlying sand bottom. Source depth is 7.3 
m, target depth is 10.3 m. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSION 

• Because of the multipath nature of acoustic 
propagation inherent in shallow water, active sonar 
acoustic pulses of finite duration are stretched in 
time leading to a significant reduction in the peak 
amplitude of the returning echo. This study 
examined the one-way energy reduction due to time 
spreading, termed energy spreading loss (ESL) , as a 
function of range, depth, source and target depth, 
bottom sediment composition and sound speed profile 
(SSP) shape. In contrast to deep water active sonar 
propagation, ESL in shallow water was found to be 
large (in excess of 10 dB for certain SSP and 
bottom sediment configurations) and exert a 
significant degradation on active sonar 
performance . 

• The time- stretch transmitted pulse propagates in a 
shallow water wave guide as discrete packets or 
modes. The shape of the time- stretched pulse is far 
from the Gaussian time distribution assumed in most 
previous analyses . 

• ESL is not a linear function of range, but 
increases rapidly out to a critical range {- 1600 m 
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for the 64 m shallow area in this study) that is 
site dependent. Beyond this critical range ESL 
remains relatively constant in the mean but 
undergoes large-scale fluctuations (standard 
deviation of 1 - 2 dB) due to modal interference 

patterns. ESL is low for those range/depth 
combinations where most of the energy is carried by 
a few dominant modes propagating with nearly 
similar group speeds. 

ESL was determined to be predominantly dependent on 
the bottom sediment composition. ESL was large for 
propagation over highly reflective (sandy) bottoms 

(8-9 dB) but moderate over absorptive (silt/clay) 

bottoms (4-5 dB) . The difference of 4 - 5 dB in 
ESL is related to the amounts of energy reflected 
from the bottom interface back into the water 
column. Highly reflective bottoms permit the 
propagation of both low and high order modes and 
their resultant variation in group speed leads to 
significant amounts of time stretching. Slow speed 
(i.e., silt/clay) absorptive bottoms attenuate the 
higher order modes (large angle rays which interact 
with the ocean boundaries) , leaving only a few low 
order modes to propagate in the water column. The 
near similar group speed of these modes results in 
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minimal to moderate time stretching. 

• In shallow water the magnitude of ESL is only- 

weakly dependent on the SSP shape, being 1 ~ 2 dB 
larger for profiles exhibiting a negative, downward 
refracting, gradient. This weak dependency on 
profile shape is a characteristic feature of 
shallow water propagation because acoustic 
interaction with the upper and lower boundaries of 
the relatively narrow wave guide will occur 
regardless of the profile shape. An implication of 
this feature is that, for any given region, ESL is 
relatively invariant with season, an important 
tactical consideration, especially for strategic 
locations where current or historical SSP 
information may be limited or lacking. 

• The transmission loss, when modeled for a single 
frequency (e.g., 3500 Hz in this study), exhibits 
rapid fluctuations of 10 ~ 20 dB along the entire 

propagation path due to phase interference of the 
propagating modes. However, when a pulse of finite 
bandwidth (200 Hz at 1 Hz increments centered on 
3.5 kHz for this study) is modeled, the 
interference pattern associated with each frequency 
tends to average or smooth the summed TL leading to 
a TL curve with virtually no spatial fluctuations. 
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Overall, a higher detection probability can be 
anticipated from a broadband source compared to 
single frequency source. 

• Modeling ESL accurately in shallow water allows one 
to quantitatively define both ESL and TL uniquely 
in order to assess the impact of the environment on 
ESL. 

• ESL plus TL or total loss must be considered in 
assessing the impact of the environment on tactical 
active sonar performance. For example, 

- ESL is low and TL is high for shallow water areas 
with silt/clay sedimentary layers. 

- ESL is high and TL is low for shallow water areas 
with hard sand bottoms. 

- considering total loss (ESL + TL) , tactical sonar 
performance is expected to be best over sand 
bottoms . 

B . RECOMMENDATIONS 

• Perform a similar analysis of the impact of the 
shallow water environment on ESL for transmitted 
pulses from advanced tactical active sonars 
(AN/SQS-53C and AN/SQS-22 (Active Low Frequency 
Sonar (ALFS) ) ) . 

• Develope advanced signal processing methods, such 
as Inverse Beam Forming (IBF) , to reduce the 
degradation incurred by ESL. 
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